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MEASUREMENT  OF  RUBIDIUM  NUMBER  DENSITY 
UNDER  OPTICALLY  THICK  CONDITIONS 


M.  ROTONDARO,  PH.D.,  LT  COL  C.  WISNIEWSKI,  AND  G.  HAGER,  PH.D. 


Abstract.  A  measurement  of  rubidium  number  density  under 
optically  thick  conditions  has  been  demonstrated  by  measuring  the 
wings  of  the  D±  absorption  spectra  using  a  laser  with  a  0.16  nm 
(75  GHz)  fine  tuning  range.  This  technique  can  measure  the  ab¬ 
solute  concentration  in  rubidium  under  conditions  where  the  ab¬ 
sorption  coefficient  and  path  length  product  yield  conditions  where 
the  central  region  of  the  line  is  opaque.  The  laser  was  tuned  to  a 
region  sufficiently  far  into  the  short  wavelength  wing  of  the  absorp¬ 
tion  where  transmission  through  the  cell  was  possible.  The  laser 
was  then  scanned  through  the  central  opaque  region  of  the  line 
to  the  adjacent  long  wavelength  wing.  The  wavelength  of  the  scan 
was  calibrated  by  using  a  1.5  GHz  etalon  and  a  cell  containing  only 
naturally  occurring  rubidium  as  a  frequency  reference.  The  mea¬ 
sured  absorption  spectra  for  various  cell  conditions  of  temperature 
and  pressure  were  then  fit  to  a  pressure  broadened  Voigt  profile 
thereby  allowing  the  determination  of  the  rubidium  number  den¬ 
sity. 


1.  Background 

In  recent  years,  alkali  metals  have  garnered  a  great  deal  of  inter¬ 
est  as  a  laser  medium.  The  small  energy  defect  between  the  D\  and 
D-2  lines  make  potassium,  rubidium  and  caesium  excellent  candidates 
for  efficient,  high-power  laser  systems.  While  these  alkali  metals  offer 
great  promise,  there  are  several  issues  which  need  to  be  resolved.  Two 
such  issues  are  heat  dissipation  and  the  challenge  of  holding  the  alkali 
number  density  constant  during  laser  operation.  This  work  will  make 
inroads  into  both  of  these  issues. 

The  first  issue  is  maintaining  a  constant  alkali  metal  concentration 
during  operation.  The  ability  to  maintain  constant  alkali  metal  concen¬ 
tration  is  dependent  on  the  ability  to  measure  the  concentration.  This 
measurement,  while  in  principle  is  not  difficult,  is  extremely  difficult 
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under  the  optically  thick  conditions  at  which  an  alkali  metal  laser  oper¬ 
ates.  The  measurement  of  the  concentration  is  performed  by  scanning 
a  laser  across  the  D\  or  D-2  line  and  measuring  the  absorption  [9].  Once 
this  absorption  spectra  is  measured  then  equation  1.1  can  be  applied 
to  relate  the  observed  spectra  to  the  number  density  N.  The  measured 
spectra  is  a  direct  measurement  of  where  a(u)  is  the  absorption 
cross-section  and  L  is  the  length  of  the  cell. 

(1.1)  In  =  —Ncr(v)L 

Under  conditions  favorable  for  the  operation  of  an  alkali  laser  the 
number  density  of  the  alkali  metal  is  sufficiently  high  to  render  the  cell 
opaque  at  line  center.  This  issue,  in  principle,  can  be  circumvented 
by  starting  a  laser  scan  in  a  wing,  scanning  through  the  opaque  region 
and  then  capturing  the  wing  on  the  other  side  of  the  peak.  These  peak 
fragments  can  then  be  fit  using  a  highly  constrained  fitting  function 
thereby  reconstructing  the  spectra.  While  this  is  possible  at  low  pres¬ 
sures  at  high  pressure  where  the  line  becomes  significantly  broad  this 
becomes  even  more  difficult.  An  alkali  laser  requires  a  high  pressure 
of  buffer  gas  to  induce  spin-orbit  relaxation  of  the  ZU  to  D\  line.  The 
buffer  gas  also  has  the  effect  of  broadening  the  lines.  Therefore,  to  ac¬ 
complish  a  wing  to  wing  laser  scan  the  laser  must  have  a  tuning  range 
on  the  order  of  75  GHz  with  a  step  size  in  the  range  of  0.5  GHz.  For¬ 
tunately,  a  laser  is  available  that  can  scan  such  a  broad  region  and  will 
be  employed  to  make  these  concentration  measurements. 

The  second  issue  is  the  heat  dissipation  resulting  from  the  spin  orbit 
relaxation  of  the  D-2  excited  state  to  the  D\  upper  laser  level.  This 
heat  deposition  along  the  laser  path  increases  the  temperature  of  the 
alkali  metal  and  buffer  gas  in  the  laser  gain  region.  This  has  the  effect  of 
degrading  the  laser  performance  [5].  To  date,  the  actual  mechanism  for 
this  performance  drop-off  has  not  been  determined.  With  the  ability  to 
measure  alkali  metal  concentration  under  conditions  favorable  for  alkali 
metal  laser  operation  it  will  be  possible  to  probe  the  gain  media  and 
make  a  determination  as  to  the  cause  of  this  performance  degradation. 

2.  Experimental  Setup 

The  experimental  apparatus  is  depicted  in  figure  1  and  a  list  of  the 
components  can  be  found  in  table  1  Starting  at  the  laser,  the  broadband 
tunable  laser  source  is  routed  through  three  beam  splitters  and  into  an 
etalon.  The  first  beam  splitter  directs  the  laser  beam  into  the  sample 
heat  pipe.  This  beam  splitter  is  situated  upon  a  translation  stage  along 
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with  it’s  paired  detector  to  facilitate  radial  dependent  number  density 
measurements  in  the  heat  pipe.  The  heat  pipe  was  enclosed  in  a  heat 
block  and  a  set  of  cold  blocks.  Heat  blocks  of  various  lengths  were 
constructed  to  facilitate  the  heating  of  different  lengths  of  the  heat  pipe. 
By  adjusting  the  heated  region  this  effectively  set  the  length  of  the  heat 
pipe.  The  heat  blocks  were  heated  using  standard  resistive  heating  and 
the  cold  blocks  were  cooled  using  a  refrigerated  water  circulator.  The 
pressure  and  composition  of  the  diluent  within  the  heat  pipe  could 
also  be  adjusted  using  the  attached  gas  handling  system.  The  gas 
handling  system  consisted  of  a  vacuum  pump,  10  Torr  and  1000  Torr 
baratrons,  various  valves  and  a  line  going  to  a  regulated  gas  cylinder. 
The  second  beam-path  was  directed  into  a  second  heat  pipe.  This 
heat  pipe  contained  only  naturally  occurring  rubidium.  This  second 
heat  pipe  was  used  as  an  absolute  wavelength  reference.  The  third 
beam-path  monitored  the  laser  power  throughout  each  scan.  This  was 
absolutely  necessary  due  to  the  large  changes  in  laser  power  during  a 
scan.  The  final  path  is  to  an  etalon.  This  etalon  was  used  to  linearize 
the  scan  axis  and  to  establish  the  correct  spacing  between  each  data 
point.  The  axis  was  then  shifted  to  the  correct  wavelength  using  the 
rubidium  wavelength  reference. 


Laser 


\ 

\  \ 

Etalon 

\ 

\ 

Sample 
Heat  Pipe 

Heat  /  Cold 
Blocks 


Translation 
stage 


Gas 

Handling 

■■ 

Rutidium 

Frequency 

Reference 


Detectors 


Laser  Power 
Reference 


FIGURE  1 .  Experimental  apparatus 
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Item 

Manufacturer 

Model 

Tunable  Laser  (794-804nm) 

New  Focus 

Velocity  6314 

Detectors  (400-1070nm) 

New  Focus 

2031 

Spectrum  Analyzer 

Coherent 

240 

Beam  Splitters 

Mclles  Griot 

BBS-485-529-1025-45 

Items  Not  Depicted 

Lock-in  Amplifiers 

Stanford  Research 

SR850 

Variable  Neutral  Density  Filter 

Mclles  Griot 

CNDQ-2-2.00 

Refrigerated  Water  Circulator 

Thermo  Scientific 

RTE7 

Vacuum  Pump 

Cole-Parmer 

79203-30 

Baratrons 

MKS 

102A,  690A 

Table  1.  Apparatus  equipment  list 


The  temperature  of  each  heat  pipe  was  monitored  using  a  thermo¬ 
couple  set  in  the  heat  block  such  that  its  sensor  was  situated  just  below 
the  bottom  of  the  heat  pipe.  This  sensor  allowed  course  control  of  the 
heat  pipe  temperature  and  allowed  monitoring  such  that  the  temper¬ 
ature  of  the  heat  pipes  could  be  held  constant  during  data  acquisition 
runs.  Although,  these  temperature  measurements  did  not  represent  a 
true  measure  of  the  temperature  of  the  alkali  metal.  This  discrepancy 
is  attributed  to  the  thermal  conductivity  of  the  heat  pipe  casing  and 
the  hot  /  cold  regions  surrounding  the  probe  region  induced  by  the 
heat  /  cold  blocks  respectively. 

The  data  was  acquired  using  two  lock-in  amplifiers.  The  two  most 
important  measurement  were  the  laser  power  and  the  sample  cell  there¬ 
fore,  these  were  acquired  using  the  phase  averaging  inputs  of  the  lock- 
ins.  The  rubidium  reference,  and  the  etalon  were  also  acquired  using 
the  lock-ins  but  these  were  recorded  using  the  auxiliary  inputs.  All 
data  was  saved  in  its  raw  format  to  allow  post  processing  of  the  data. 

The  laser  has  a  large  tuning  region  (794-804  nm)  which  can  be 
scanned,  using  the  stepper  motor,  in  increments  as  small  as  0.01  nm 
(4.75  GHz).  This  step  size  is  too  coarse  for  the  purpose  of  measuring  the 
pressure  broadening  and  therefore  the  number  density.  Subsequently, 
it  was  necessary  to  use  the  fine  frequency  scan  mode  of  the  laser  which 
employs  a  piezoelectric  crystal.  The  piezoelectric  crystal  position  as 
a  function  of  input  voltage  can  be  seen  if  figure  2.  Additionally,  the 
discrete  nature  of  the  piezoelectric  crystals  steps  can  be  seen  in  the 
blown-up  region  on  the  figure.  The  piezoelectric  crystal  voltage  could 
be  set  between  3.8%  and  99.8%  of  maximum  in  increments  of  0.4%. 
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The  fine  region  scan  length  is  0.16  nm  (75G  Hz)  therefore,  the  min¬ 
imum  fine  scan  step  size  is  approximately  0.67  pm  (315  MHz).  This 
scan  size  and  the  finesse  of  the  etalon  put  a  limit  on  the  highest  reso¬ 
lution  etalon  that  could  be  used  because  of  the  necessity  to  access  all 
etalon  resonance  peaks.  Therefore,  a  (1.5  GHz)  etalon  was  chosen  as  a 
wavelength  reference. 


Input  Piezoelectric  Crystal  Voltage 


FIGURE  2.  Piezoelectric  crystal  position  as  a  function  of 
the  input  voltage  to  the  crystal.  The  emphasized  region 
indicates  the  discrete  nature  of  the  piezoelectric  crystals 
steps. 


3.  Data  Acquisition 

In  preparation  for  taking  data,  the  first  step  was  to  ensure  the  heat 
pipes  were  set  to  the  appropriate  temperature.  For  the  sample  heat 
pipe,  this  meant  waiting  a  sufficient  length  of  time  for  the  temperature 
of  the  heat  block  to  fully  stabilize.  For  the  reference  cell,  this  meant 
ensuring  there  was  a  sufficiently  strong  absorption  spectra  to  ensure  a 
good  quality  fit  later  in  the  analysis. 

During  the  warmup  period  for  the  heat  pipes,  the  laser  was  set  to  its 
initial  scanning  position.  This  was  accomplished  by  setting  the  initial 
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wavelength  of  the  laser  to  a  sufficiently  low  wavelength  to  ensure  it  was 
well  outside  of  any  of  the  absorption  features.  Then  the  piezoelectric 
crystal  was  set  to  50%  of  its  range.  Then,  the  laser  was  scanned  slowly, 
using  the  laser’s  stepper  motor  on  its  smallest  step  size,  towards  longer 
wavelengths  until  the  reference  cell  absorption  spectra  began  to  appear 
on  the  data  acquisition  monitor.  At  that  point  the  laser’s  course  scan 
was  stopped.  This  placed  the  laser  at  the  center  of  the  absorption  fea¬ 
ture.  The  laser  could  then  be  scanned  from  0-100%  of  the  piezoelectric 
crystal  range  allowing  the  acquisition  of  the  absorption  spectra.  The 
laser  can’t  really  scan  from  0-100%  it  is  actually  limited  to  3.8-99.8%. 
By  inverting  the  data  represented  in  figure  2  and  fitting  the  result  to  a 
line.  It  is  possible  to  improve  the  linearity  of  the  scan.  The  resulting 
fit  is  given  in  equation  3.1. 


(3.1)  VA ct  =  1.03067Vbe*  -  3.013635 

Where  Vdct  is  the  actual  voltage  sent  to  the  piezoelectric  crystal  and 
Vfjes  is  the  desired  voltage  setting.  Vjjes  can  range  from  3.8-99.8%.  This 
correction  to  the  laser  scan  is  not  really  necessary  because  any  nonlin¬ 
earity  in  the  scan  will  be  corrected  using  the  data  from  the  etalon.  The 
reason  for  the  application  of  this  correction  is  to  reduce  the  difficulty 
of  fitting  the  etalon  data. 

When  the  laser  is  scanned  four  channels  of  data  are  collected.  These 
channels  are  the  signal  from  the  sample  cell,  rubidium  reference  cell, 
laser  power  and  the  etalon.  The  data  was  acquired  using  two  lock-in 
amplifiers.  The  two  most  important  measurement  were  the  laser  power 
and  the  sample  cell.  The  laser  power  because  all  other  measurements 
were  ratioed  against  it  to  correct  for  the  laser’s  power  fluctuations 
during  the  scan.  The  sample  cell  because  it  is  the  data  used  to  measure 
the  rubidium  number  density.  Therefore,  these  were  acquired  using 
the  phase  averaging  inputs  of  the  lock-in  amplifiers.  The  rubidium 
reference,  and  the  etalon  were  also  acquired  using  the  lock-in  amplifiers 
but  these  were  recorded  using  the  auxiliary  inputs.  All  data  was  saved 
in  its  raw  format  to  allow  post  processing  of  the  data.  Typical  raw  data 
can  be  seen  in  figure  3.  As  can  be  seen,  the  laser  power  has  a  large 
variation  over  the  entire  scan  region.  This  fluctuation  is  corrected  by 
ratioing  all  of  the  collected  data  by  the  laser  power.  The  ratioed  data 
can  be  seen  in  figure  4.  This  data  represents  —  ln(/sa™pie ),  —  ln(-^) 
and  lEt^on .  As  can  be  seen  the  fluctuations  in  the  laser  power  are 
corrected. 


7 


1.4 
1.2  - 
1.0  - 
0.8  - 
0.6  - 
0.4  - 
0.2  - 
0.0  - 
-0.2 
-0.4 
-0.6 
-0.8 


Data  Cell 
Rb  Reference 
Laser  Power 


-60  -40  -20  0  20 

Oscilloscope  Time  Step 


40 


60 


Figure  3.  Typical  raw  data 


The  X-axis  of  the  collected  data  represents  the  divisions  used  by 
the  data  acquisition  equipment  and  needs  to  be  converted  into  a  wave¬ 
length.  This  is  accomplished  by  using  the  etalon  data.  The  procedure 
for  using  the  etalong  data  is  explained  it  detail  in  [7]  but  to  summarize. 
The  etalon  data  is  fit  using  the  software  package  PeakFit  Version  4.12. 
This  software  is  very  good  at  finding  the  peak  center  of  a  long  line  of 
successive  peaks.  It  is  critical  that  the  software  identify  every  peak. 
Once  each  peak  is  identified,  an  X-Y  table  is  created  using  the  peak 
centers  as  the  X  values  and  the  incremental  spacing  of  each  peak  as 
determined  by  the  etalon  spacing  for  the  Y  values.  Once  this  table  is 
created  it  is  then  fit  to  a  third  degree  polynomial.  This  polynomial  is 
then  used  to  correct  for  any  nonlinearity  and  to  set  the  correct  spacing 
between  adjacent  points.  This  results  in  an  x-axis  that  has  the  correct 
wavelength  spacing  but  has  not  been  shifted  to  the  correct  wavelength. 
This  correction  is  applied  to  both  the  sample  cell  data  and  the  rubidium 
reference  cell  data. 

The  next  step  is  to  shift  the  x-axis  to  the  correct  wavelength.  This 
is  accomplished  by  fitting  the  rubidium  reference  data  to  a  series  of 
Gaussian  curves.  The  gaussian  line  shape  is  chosen  because  the  rubid¬ 
ium  is  dominated  by  Doppler  broadening  which  results  in  a  Gaussian 
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Figure  4.  Data  that  has  been  divided  by  the  laser  power 


line  shape.  The  relative  spacing  and  line  strengths  of  the  rubidium 
lines  has  been  well  documented  [10]  [11]  and  the  exact  Di  line  shape 
g(X)  or  g{y)  can  be  reconstructed  using  equations  3.2  or  3.3 

(3.2)  g(A)  =  f F  f isoS FF'  *  9g  XfF') 

F,F'  , iso 


(3.3) 


9(0=  J2  f Ffiso^FF’  *  QGaussiy i  VFF>) 

F,F',iso 

Where  the  Gaussian  in  both  wavelength  and  frequency  are  given  by 
equations  3.4  and  3.5. 

,  .  x  ,  /41n(2)  1  r  /a  a  x  21 

(3.4)  gGauss{X,  Aff'J  =  V — I — WR- exP 


7 r  AA 


D 


,  .  ,  .  /41n(2)  1 

(3.5)  gcaussiv,  VFF>)  =  \j  — - — -x—  exp 


7 r  Av 


D 


-4  ln(2)  | 

(A  —  Aff* 

V  AA  d 

4 ln(2)  ( 

^  v  —  Vpp> 

\  Avd 

The  two  forms  of  the  line  shape  are  related  by  equation  3.6 


(3.6) 


A2 

gGaussi A,  \FF>)  =  g Gauss ( F  Fppi) 

C 
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and  both  are  normalized  when  integrated  over  all  frequencies  as  shown 
by  equations  3.7. 

poo  POO 

(3.7)  /  gGaussi^,  \FF')d\  =  /  gGauss(v,  VFF')dlS  =  l 

Jo  Jo 

The  full  width  at  half  maximum  (FWHM)  for  each  form  is  given  in 
equation  3.8. 


(3.8) 


AAn  —  Ar 


8  kT  ln(2) 


,  A  uD  =  f0 


8 kT  ln(2) 


Me2  ’  u  V  Me2 

Where  k  is  the  Boltzmann  constant,  T  is  the  temperature,  M  is  the 
mass  of  a  rubidium  atom  and  c  is  the  speed  of  light. 

In  the  previous  equations,  F  is  the  hyperfme  quantum  number  of  the 
starting  level  and  F’  is  the  hyperfme  quantum  number  of  the  ending 
level,  fp  is  the  statistical  distribution  of  the  F  state  and  is  given  by 
equation  3.9. 


(2F  +  l)r^'g  _  (2F+1) 

1  '  lF  EA2F  +  l)e-E<F>/tT  ~  Ef(2F  +  1) 

This  approximation  is  allowed  because  the  energy  E(F)  is  very  small 
representing  the  energy  defect  between  the  hyperfme  split  levels.  Ex¬ 
amining  the  worst  case,  the  hyperfme  splitting  of  the  ground  state, 
exp(-^P-)  ~  0.999.  The  isotopic  abundance  is  represented  by  fiSO  and 
is  0.2783  for  Rubidium  87  and  0.7217  for  Rubidium  85.  The  product 
Sff'  *  If  is  given  in  table  3.  Spp>  is  the  relative  hyperfme  transition 
strength  factors  given  in  table  2.  A pp>  is  the  line  center  of  the  transition 
F  to  F’. 


multiplied  by  the  statistical  distribution  of  the  hyperhne 
levels.  Sff'  *  /f 
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The  line  centers  A ff>  are  determined  by  calculating  the  line  positions 
relative  to  the  Si(F=2)  to  Pi(F=2)  transition  in  87 Rb  and  S'i(F=3) 
to  Pi(F=2)  transition  in  85  Rb  for  the  D\  line. 

The  total  energy  of  a  hyperhne  split  level  with  the  total  angular 
momentum  F  is  given  by  equation  3.10  from  reference  [6]. 


(3.10) 

with 


C 

Ep  =  Ej  +  A— 


|C(C  +  l)-/(/  +  l)J(J  +  l) 
2/(2/  —  1)  J(2J  —  1) 


C  =  F(F  +  1)  —  /(/  +  1)  —  </(</  +  1); 

F  =  /  +  J,i  +  J-i,...|/- J\; 

_  /Q//(0), 

IJ  ’ 

B  =  eQ4>jj(  0). 


Pj  is  the  fine  structure  energy,  and  the  hyperhne  splitting  is  given 
by  the  last  two  terms.  The  term  Aj  represents  the  splitting  due  to 
the  nuclear  magnetic  moment  which  is  dependent  upon  the  magnetic 
moment  Hi,  the  magnetic  held  produced  by  the  nucleus  H( 0)  and  the 
product  of  the  quantum  numbers  IJ  which  represents  the  coupling 
of  the  angular  momentum  vectors  through  the  magnetic  interaction 
where  I  —  |  for  8‘  Rb  and  I  —  |  for  85 Rb.  The  last  term  represents  the 
quadrupole  interaction  which  depends  on  Q,  the  quadrupole  moment, 
(f)jj  (0) ,  the  vector  gradient  of  the  electric  held  of  the  orbital  electrons 
having  cylindrical  symmetry  about  the  J  axis  evaluated  at  the  nucleus, 
and  e  is  the  absolute  value  of  the  elementary  charge. 

The  relative  line  positions  are  given  in  Table  4.  The  hrst  half  of 
the  table  provides  the  A  and  B  spectroscopic  constants  which  can  be 
used  to  determine  the  relative  line  positions  as  in  Figure  5.  The  figure 
accurately  represents  the  relative  line  positions  within  each  manifold 
but  the  spacing  between  manifolds  is  not  to  scale.  The  relative  line  po¬ 
sitions  are  computed  in  the  second  half  of  Table  4  using  equation  3.11. 


(3.11) 


rpFF' 


All  line  positions  are  reported  relative  to  the  Si  (F  =  2)  to  Pi  (. F  =  2) 
transition  in  87 Rb  for  the  Di  line  and  to  the  Si(F  =  2)  to  Ps(F  =  3) 
transition  in  81  Rb  for  the  D2  line.  The  85 Rb  shifts  are  relative  to  the 
Si  (F  =  3)  to  Pi  (F  =  2)  transition  for  the  Di  line  and  to  the  Si  (. F  =  3) 
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87  Rb  Ref  [1]  [2]  [3]  [4]  [12] 

85  Rb  Ref  [1]  [2]  [3] 

MHz 

A 

B 

A 

B 

2 

3417.341305452145(45) 

- 

1011.9108130(20) 

- 

5  2Pi 

2 

407.24(77) 

- 

120.527(56) 

- 

5  2Pl 

2 

84.7185(20) 

12.4965(37) 

25.0020(99) 

25.790(93) 

rpi 

87  Rb 

85  Rb 

nm 

F 

1 

2 

F 

2 

3 

Pi 

2 

1 

-0.0126821 

0.0017159 

2 

-0.0063951 

0 

2 

-0.014398 

0 

3 

-0.0071568 

-0.0007617 

P3 

2 

0 

-0.013353 

- 

1 

-0.006333 

- 

1 

-0.013506 

0.000892 

2 

-0.006395 

0 

2 

-0.013836 

0.000562 

3 

-0.006528 

-0.0001335 

3 

- 

0 

4 

- 

-0.0003877 

GHz 

F 

1 

2 

F 

2 

3 

Pi 

2 

1 

6.02014 

-0.81454 

2 

3.03573 

0 

2 

6.83468 

0 

3 

3.39731 

0.361581 

Ps 

2 

0 

6.33886 

- 

1 

3.00636 

- 

1 

6.41109 

-0.423592 

2 

3.03573 

0 

2 

6.56803 

-0.266652 

3 

3.09913 

0.06340 

3 

- 

0 

4 

- 

0.18404 

cm  1 

F 

1 

2 

F 

2 

3 

Pi 

2 

1 

0.20081 

-0.02717 

2 

0.10126 

0 

2 

0.22798 

0 

3 

0.11332 

0.01206 

Ps 

2 

0 

0.21144 

- 

1 

0.10028 

- 

1 

0.21385 

-0.014129 

2 

0.10126 

0 

2 

0.21908 

-0.008894 

3 

0.10337 

0.00211 

3 

- 

0 

4 

- 

0.00614 

Table  4.  A  and  B  coefficients  and  relative  line  posi¬ 


tions,  rpi 


to  P3(F  =  2)  transition  for  the  D2  line.  This  requires  an  additional  pa¬ 
rameter,  the  spacing  between  the  81  Rb  and  the  85 Rb  lines.  This  spacing 
has  been  calculated  using  the  values  reported  in  [10]  [11]  for  the  Di  line 
as  —0.0014843  nm  (704.07  MHz)  (0.023485  cm-1)  and  for  the  D2  line 
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as  —0.0019138  nm  (942.44  MHz)  (0.031436  cm-1).  The  absolute  wave¬ 
length  correction  of  the  axis  requires  that  these  transitions  be  shifted 
to  the  absolute  position  of  the  Si(F  =  2)  to  Pi(F  =  2)  transition  in 
87 Rb  which  is  794.9836  nm  (377.1059  THz)  and  for  the  Si(F  =  2)  to 
Ps(F  =  3)  transition  in  87 Rb  it  is  780.246  nm  (384.228  THz).  Once 
these  values  have  been  calculated  the  line  shape  can  be  assembled  using 
equation  3.2  and  this  can  be  seen  graphically  in  figure  6. 

In  order  to  extract  the  shift  needed  to  correct  the  data’s  X-axis,  this 
constructed  line  shape  is  assembled  as  a  user  function  in  the  software 
package  TableCurve  2D  v5.01.01.  The  etalon  corrected  reference  data 
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794.96  794.97  794.98  794.99  795.00 

Wavelength  (nm) 


FIGURE  6.  Natural  rubidium  with  each  D\  hyperfine 
transition  explicitly  shown 

is  then  fit  allowing  the  extraction  of  the  shift  necessary  to  correct  the 
wavelength  axis  and  the  Doppler  width.  Reference  data  with  the  its 
resulting  curve  fit  are  shown  in  figure  7.  The  Doppler  width  is  used  to 
determine  the  temperature  of  the  rubidium  within  the  heat  pipe. 

By  using  the  fit  of  the  etalon  peaks  and  the  fit  from  the  rubidium 
reference,  the  data  can  be  both  scaled  and  shifted  such  that  the  x-axis 
represents  the  wavelength  of  the  spectra.  There  is  an  additional  step 
needed  to  correct  the  axis.  The  etalon  is  actually  a  spectrum  analyzer 
therefore,  one  of  it’s  mirrors  is  attached  to  a  piezoelectric  crystal.  This 
crystal  can  move  from  the  application  of  an  external  voltage  or  from 
changes  in  room  temperature.  Therefore,  it  is  necessary  to  do  a  second 
set  of  fits  on  the  data.  The  etalon  peaks  are  fit  and  the  resulting 
fit  is  applied  to  the  reference  data.  The  reference  data  is  then  fit 
allowing  the  x-axis  to  scale.  The  scale  parameter,  is  then  applied  to  the 
spacing  of  the  etalon  peaks  and  the  etalon  peak  spacing  is  readjusted 
to  account  for  this  small  correction.  Then,  the  result  is  applied  to  both 
the  reference  cell  and  the  sample  cell  data.  The  resulting  reference  data 
is  then  fit  one  final  time  to  ensure  the  x-axis  has  been  properly  scaled 
and  shifted.  Data  with  the  wavelength  axis  corrected  is  depicted  in 
figure  8 


14  M.  ROTONDARO,  PH.D.,  LT  COL  C.  WISNIEWSKI,  AND  G.  HAGER,  PH.D. 


Figure  7.  Rubidium  87  data  with  fits 


4.  Data  Analysis 


Once  the  data  has  been  acquired  and  the  wavelength  axis  has  been 
scaled  and  shifted  to  the  correct  wavelength  region  it  can  be  used  to 
extract  the  rubidium  number  density.  The  sample  heat  pipe  contained 
naturally  occurring  rubidium  and  was  Filed  to  various  pressures  ranging 
between  0  Torr  and  1000  Torr  using  helium,  nitrogen  and  ethane  as 
diluent.  Because  of  the  addition  of  the  diluent  and  the  various  pressure 
ranges,  the  spectra  in  the  sample  cells  consisted  of  both  a  Gaussian  line 
shape  and  a  pressure  broadened  Lorentzian  line  shape.  Therefore,  the 
resulting  line  shape  is  composed  of  the  sum  eight  Voigt  profiles.  A 
Voigt  line  shape  is  represented  by  equations  4.1  and  4.2. 

(4.1) 


9  V  oigt  (  A  5  A  FF'J 


AXl  exp (— t2) 

<A  -  V7  -  +W 


dt 


1 

r  oo 

A uL  exp(— t2) 

2n^  j 

L 

V  —  ^)2+(V)2J 

(4.2) 

9Voigt{v ,  UfF') 


dt 
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Wavelength  (nm) 


Figure  8.  Data  that  has  been  scaled  and  shifted  using 
the  etalon  and  reference  data 

Using  the  Voigt  line  shape  and  equation  3.2  the  line  shape  for  the 
rubidium  sample  heat  pipe  spectra  can  be  constructed.  The  resulting 
Voigt  line  shape  can  be  curve  fit  to  the  data  and  the  Voigt  full  width 
at  half  maximum  (FWHM)  and  the  amplitude  can  be  obtained.  This 
method  works  well  when  the  entire  rubidium  spectra  is  observed  but 
under  highly  opaque  conditions  the  Voigt  profile  is  quite  difficult  to 
fit.  To  improve  the  quality  of  the  fit  and  subsequently  the  quality  of 
the  number  density  measurement,  the  Gaussian  and  the  Lorentzian 
FWHM  are  calculated  using  equations  3.8  and  4.3  respectively. 

(4.3)  AAl  =  R(X)P,  A vL  =  R{y)P 

Where  P  is  the  pressure  of  the  diluent  and  R  is  the  broadening  rate 
given  in  or  The  Ethane  value  used  for  the  broadening  rate  R 

was  27.8  measured  at  314.15  K  and  it  was  obtained  from  [13].  The 

Helium  value  used  for  the  broadening  rate  R  was  18.90  yp^  measured 
at  394  K  and  it  was  obtained  from  [8]. 

The  Doppler  FWHM  is  temperature  dependent  and  therefore  can’t 
be  calculated  without  a  precise  measurement  of  the  temperature.  That 
temperature  can  be  measured  by  measuring  the  Doppler  width  of  the 
rubidium  at  0  Torr  of  diluent.  While  the  temperature  may  change  as 
diluent  is  added  to  the  heat  pipe  the  y/T  dependance  of  the  doppler 
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width  will  vary  slowly  and  therefore  can  be  neglected.  These  FWHM 
values  are  then  fixed  and  the  only  parameters  that  are  allowed  to  float 
are  the  baseline  and  the  amplitude  of  the  Voigt  profile.  Typical  data 
with  their  associated  fits  is  depicted  in  figure  9. 
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Figure  9.  Rubidium  87  data  with  fits 


Once  the  amplitude  of  the  spectra  is  determined  Beer’s  Law  equa¬ 
tion  1.1  and  the  relationship  between  the  absorption  cross  section 
and  the  line  shape  equation  4.4  can  be  combined  to  for  equations  4.9 
and  4.10  which  are  used  to  calculate  the  number  density.  This  analysis 
was  conducted  using  wavelength  as  the  X-axis  therefore,  the  wavelength 
form  of  the  equations  were  used.  The  frequency  form  of  the  equations 
were  included  for  completeness.  The  results  of  these  calculations  can  be 
seen  in  figure  10  also  depicted  is  the  number  density  calculation  using 
the  vapor  pressure  curve  equation  4.11  from  [10]  which  is  reported  to 
have  an  accuracy  better  than  ±5%.  Note,  this  equation  has  been  mod¬ 
ified  to  yield  number  density  N  /m3  as  a  function  and  temperature  T 
in  Kelvin. 
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(4.4) 

A2  A4 

cr{y)  =  - — g(i/,  uFF/,  AvD,  AuL),  cr(A)  =  - - g(X,  XFF>,  AXd,  AXl) 

o7TT  87 TTC 


(4.5)  -  In  (  =  Na(v)L  =  NL^^g(is,  vFF>,  AuD,  AvL ) 

h\y) )  ottt 


(4.6)  -  In  (  =  Na(X)L  =  NL-^—g(X,  \FF,,  AXD,  AXL) 

loyA)  J  OTTTC 


A2 

(4.7)  AFit(u)g(u,  vfffAvd,  A uL)  =  NL- — g(v,  vFF>,  AvD ,  A uL) 

OTTT 

A4 

(4.8)  AFit(X)g(X,  XFF>,  AXF,  AX F)  =  NL- - <7 (A,  XFF>,  AXF,  AX F) 

87 TTC 


(4.9) 

(4.10) 

(4.11) 


N  —  AFit{y) 
N  =  Aplt(X) 


8ttt 

LP 

87 TTC 

Tx* 


2.077  x  109  ^ 

N  = - — - •  10  t 

kT 


Note,  Apitiy)  is  in  units  of  1/time  and  AFu{ A)  is  in  units  of  length. 


5.  Conclusions 

As  can  be  seen  from  figure  10  there  is  very  good  agreement  between 
the  number  density  measured  by  directly  probing  the  heat  pipe  when 
compared  to  the  expected  value  derived  from  equation  4.11  using  the 
temperature  measured  from  the  doppler  line  shape.  The  measurements 
at  higher  temperatures  resulted  in  a  larger  temperature  uncertainty 
because  the  rubidium  vapor  pressure  was  sufficiently  high  to  render 
the  cell  opaque  at  line  center  of  the  hyperffiie  split  spectra.  There¬ 
fore,  the  0  Torr  buffer  gas  spectra  was  fit  using  peak  fragments  as  the 
high  temperature  number  densities  are  measured.  This  resulted  in  a 
greater  temperature  uncertainty.  Even  with  this  uncertainty,  the  num¬ 
ber  density  measurement  is  accurate  because  of  the  weak  temperature 
dependance  and  small  contribution  to  the  line  shape  of  the  Doppler 
portion  of  the  Voigt. 
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Figure  10.  Measured  rubidium  number  density  as  a 
function  of  temperature  along  with  the  vapor  pressure 
curve  for  rubidium. 


This  methodology  demonstrates  a  robust  capability  for  measuring 
the  number  density  of  rubidium  under  conditions  that  are  favorable 
for  rubidium  laser  operation.  Therefore,  the  application  of  this  tech¬ 
nique  to  diode  pumped  alkali  metal  lasers  (DPAL)  would  yield  valuable 
information  that  would  enhance  both  our  theoretical  and  experimental 
knowledge  of  DPAL  systems. 


6.  Appendix  1 

This  is  the  curve  fitting  formalism  used  to  perform  the  curve  fits. 
This  formalism  develops  the  curve  fit  using  an  approximation  of  the 
Voigt  profile.  The  actual  curve  fitting  process  used  equation  6.1  which 
was  integrated  using  an  adaptive  Gauss  Quadrature. 


1 

poo 

Az/Lexp(— t2) 

2ny/n  j 

—  OO 

(u  vFFt  vs  t  )2  +  (AD2 

_V  S  y/  4  ln(2) '  V  2  '  _ 

(6.1) 

V (u,  VFF')  =  OL 


dt+b 
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(6.2)  G(v,vFF')  =  a\/  4Ill(2)  -J—  exp 


(6.3) 

(6.4) 

(6.5) 


7T  Az/p 
Zy(z/,  UpF')  —  01 


— 41n(2) 
A  vL 


V  —  Vppi  —  Vs 


Au 


D 


^[(^f^)2  +  ("  -  VpF’  ~  Vs)2} 


U(u,  vFF')  =  (1  -  r])G(v,  vFF')  +  r)L(v ,  z/FjF/)  +  6 


V  = 


a 2  + 1 


/g  g-,  _  Avp  i  _ _ 2 _ 

AV  1  +  0.0990  ln(2)  +  [(1  -  0.0990  ln(2))2  + 
Nonlinear  least  squares  fit 

(6.7)  dAij  = 


(6.8) 

(6.9) 

(6.10) 


dV 

l^l 

dV 

l^l 

dU 

dU 

dU 

1 

da 

db 

dAuL 

Wi 

dAvn 

\V\ 

dva 

W  i 

8V 

\v2 

dV 

\v2 

au 

1 

dir 

1 

dU 

1 

da 

db 

9AuL 

\V2 

9Aud 

\V2 

dvs 

k2 

av  | 

dV\ 

dU  | 

dU  | 

dU  | 

da  1 

Vm 

db  1 

Vm 

9Aul  1 

Vm 

dAujy  1 

Vm 

dva  1 

(6.11) 


dBi  =  yi-V (i/;  a,  b,  AuL,  AuD,  vs) 

dATdB  =  (dATdA)dX 


A  =  dArdA.  B  =  dATdB 


da 
db 

dXi  =  |  dAvp 
dAvp, 
dp* 


Solved  using  the  Jacobi  Method 

B<  -  Ajdxt" 


(6.12) 


dx)k)  = 


A. 


l,Z 


Iterate  until  done... 
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(6.13) 

(6.14) 

(6.15) 

(6.16) 

dU 
dA  uD 

(6.17) 

dU  _ 


dV 

da 


1 

poo 

Avl  exp(— t2) 

2vt^/7t  j 

— oo 

[("  -  vff'  tjSk>)2 + (T,2J 

dV_ 

~db 
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Figure  11.  Need  Text  Here 


21 


A 

b 

A  vL 

A  vD 

vs 

1 

0.2 

1 

0.333 

0.2 

1.005 

0.201 

0.992 

0.367 

0.196 

1 

0.5 

1 

0.5 

0.5 

0.996 

0.501 

1.005 

0.487 

0.504 

1 

1 

1 

1 

1 

0.768 

1.016 

0.063 

1.646 

1.019 

1 

2 

1 

2 

2 

1.087 

1.995 

1.477 

1.659 

1.993 

1 

0.2 

0.2 

1 

-0.2 

1.138 

0.189 

0.519 

0.793 

-0.197 

1 

.5 

.5 

1 

-.5 

0.861 

0.512 

0.083 

1.267 

-0.497 

1 

1 

1 

1 

-1 

1.026 

0.997 

1.069 

0.940 

-1.015 

1 

2 

2 

1 

-2 

0.960 

2.003 

1.867 

1.148 

-2.002 

Table  5.  Comparison  of  Voigt  parameters  versus  the 
values  obtained  from  a  fit. 
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